Structural basis for the mechanism of interaction of SARS-CoV-2 B.1.640.2 variant RBD with the host receptors hACE2 and GRP78.

The inflicted chaos instigated by the SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) globally continues with the emergence of novel variants. The current global outbreak is aggravated by the manifestation of novel variants, which affect the effectiveness of the vaccine, attachment with hACE2 (human Angiotensin-converting enzyme 2) and immune evasion. Recently, a new variant named University Hospital Institute (IHU) (B.1.640.2) was reported in France in November 2021 and is spreading globally affecting public healthcare. The B.1.640.2 SARS-CoV-2 strain revealed 14 mutations and 9 deletions in spike protein. Thus, it is important to understand how these variations in the spike protein impact the communication with the host. A protein coupling approach along with molecular simulation protocols was used to interpret the variation in the binding of the wild type (WT) and B.1.640.2 variant with hACE2 and Glucose-regulating protein 78 (GRP78) receptors. The initial docking scores revealed a stronger binding of the B.1.640.2-RBD with both the hACE2 and GRP78. To further understand the crucial dynamic changes, we looked at the structural and dynamic characteristics and also explored the variations in the bonding networks between the WT and B.1.640.2-RBD (receptor-binding domain) in association with hACE2 and GRP78, respectively. Our findings revealed that the variant complex demonstrated distinct dynamic properties in contrast to the wild type due to the acquired mutations. Finally, to provide conclusive evidence on the higher binding by the B.1.640.2 variant the TBE was computed for each complex. For the WT with hACE2 the TBE was quantified to be-61.38 ± 0.96 kcal/mol and for B.1.640.2 variant the TBE was estimated to be -70.47 ± 1.00 kcal/mol. For the WT-RBD-GRP78 the TBE -was computed to be 32.32 ± 0.56 kcal/mol and for the B.1.640.2-RBD a TBE of -50.39 ± 0.88 kcal/mol was reported. This show that these mutations are the basis for higher binding and infectivity produced by B.1.640.2 variant and can be targeted for drug designing against it.Communicated by Ramaswamy H. Sarma.

Protein degradation: a novel computational approach to design protein degrader probes for main protease of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has afflicted many lives and led to approvals of drugs and vaccines for emergency use. Even though vaccines have emerged, the high mortality of COVID-19 and its insurgent proliferation throughout the masses commands an innovative therapeutic proposition for the treatment. Targeted protein degradation has been applied to various disease domains and we propose that it could be incredibly beneficial to tackle the current pandemic. In this study, we have attempted to furnish insights on the design of suitable PROTACs for the main protease (Mpro) of SARS-CoV-2, a protein that is considered to be an essential target for viral replication. We have employed protein-protein docking to predict the possible complementarity between a cereblon E3 ligase and Mpro of SARS-CoV-2, and estimate possible linker length. Molecular Dynamic simulation and analysis on generated ternary complexes demonstrated stable interactions that suggested that designed PROTAC has a potential to cause degradation. The superior characteristics rendered by PROTACS led us to propose them as possibly the next-generation antiviral drugs for SARS-CoV-2.

Understanding the molecular interaction of SARS-CoV-2 spike mutants with ACE2 (angiotensin converting enzyme 2).

Covid-19 is a viral disease caused by the virus SARS-CoV-2 that spread worldwide and caused more than 4.3 million deaths. Moreover, SARS-CoV-2 still continues to evolve, and specifically the E484K, N501Y, and South Africa triple (K417N + E484K + N501Y) spike protein mutants remain as the 'escape' phenotypes. The aim of this study was to compare the interaction between the receptor binding domain (RBD) of the E484K, N501Y and South Africa triple spike variants and ACE2 with the interaction between wild-type spike RBD-ACE2 and to show whether the obtained binding affinities and conformations corraborate clinical findings. The structures of the RBDs of the E484K, N501Y and South Africa triple variants were generated with DS Studio v16 and energetically minimized using the CHARMM22 force field. Protein-protein dockings were performed in the HADDOCK server and the obtained wild-type and mutant spike-ACE2 complexes were submitted to 200-ns molecular dynamics simulations with subsequent free energy calculations using GROMACS. Based on docking binding affinities and free energy calculations the E484K, N501Y and triple mutant variants were found to interact stronger with the ACE2 than the wild-type spike. Interestingly, molecular dynamics and MM-PBSA results showed that E484K and spike triple mutant complexes were more stable than the N501Y one. Moreover, the E484K and South Africa triple mutants triggered greater conformational changes in the spike glycoprotein than N501Y. The E484K variant alone, or the combination of K417N + E484K + N501Y mutations induce significant conformational transitions in the spike glycoprotein, while increasing the spike-ACE2 binding affinity.Communicated by Ramaswamy H. Sarma.

In Silico Analysis of Bacteriocins from Lactic Acid Bacteria Against SARS-CoV-2.

The COVID-19 pandemic caused by a novel coronavirus (SARS-CoV-2) is a serious health concern in the twenty-first century for scientists, health workers, and all humans. The absence of specific biotherapeutics requires new strategies to prevent the spread and prophylaxis of the novel virus and its variants. The SARS-CoV-2 virus shows pathogenesis by entering the host cells via spike protein and Angiotensin-Converting Enzyme 2 receptor protein. Thus, the present study aims to compute the binding energies between a wide range of bacteriocins with receptor-binding domain (RBD) on spike proteins of wild type (WT) and beta variant (lineage B.1.351). Molecular docking analyses were performed to evaluate binding energies. Upon achieving the best bio-peptides with the highest docking scores, further molecular dynamics (MD) simulations were performed to validate the structure and interaction stability. Protein-protein docking of the chosen 22 biopeptides with WT-RBD showed docking scores lower than -7.9 kcal/mol. Pediocin PA-1 and salivaricin P showed the lowest (best) docking scores of - 12 kcal/mol. Pediocin PA-1, salivaricin B, and salivaricin P showed a remarkable increase in the double mutant's predicted binding affinity with -13.8 kcal/mol, -13.0 kcal/mol, and -12.5 kcal/mol, respectively. Also, a better predicted binding affinity of pediocin PA-1 and salivaricin B against triple mutant was observed compared to the WT. Thus, pediocin PA-1 binds stronger to mutants of the RBD, particularly to double and triple mutants. Salivaricin B showed a better predicted binding affinity towards triple mutant compared to WT, showing that it might be another bacteriocin with potential activity against the SARS-CoV-2 beta variant. Overall, pediocin PA-1, salivaricin P, and salivaricin B are the most promising candidates for inhibiting SARS-CoV-2 (including lineage B.1.351) entrance into the human cells. These bacteriocins derived from lactic acid bacteria hold promising potential for paving an alternative way for treatment and prophylaxis of WT and beta variants.

South African (501Y.V2) and the United Kingdom (B.1.1.7) SARS-CoV-2 spike (S) protein variants demonstrate a higher binding affinity to ACE2.

INTRODUCTION: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), affects the lower respiratory tract by binding to angiotensin-converting enzyme 2 (ACE2) via its S-protein. Recent emerging SARS-CoV-2 variants from the United Kingdom (B.1.1.7) and South Africa (501Y.V2) are spreading worldwide at an alarming rate. The new variants have manifested amino acid substitution K417N, E484K and N501Y on the RBD domain that binds to ACE2. As such, these mutations may influence the binding of the S-protein to ACE2 and affect viral entry into the host cell. Methods In this study, we modelled the amino acids substitutions on the S-protein and utilised HADDOCK server to assess the S-protein RBD domain binding with ACE2. Additionally, we calculated the binding affinity of ACE2 to S-protein WT, B.1.1.7 and 501Y.V2 variants using Molecular Mechanics-Generalized Born Surface Area (MM/GBSA). Results We demonstrate that the S-protein of both variants possesses higher binding affinity to ACE2 than WT, with the South African 501Y.V2 is a more infective strain than the B.1.1.7 that originated in the United Kingdom. Conclusion The South African 501Y.V2 variant presents three amino acid substitutions that changed the H-bonding network resulting in a higher affinity to ACE2, indicating that the 501Y.V2 strain is more infective than the B.1.1.7 strain.

A Non-Structural Protein of Severe Acute Respiratory Syndrome Coronavirus 2 with a Potential Ability to Reduce Blood Glucose for Use in Controlling Type 2 Diabetes Mellitus.

It is estimated that more than 400 million people worldwide are suffering from diabetes. There are two types of diabetes. Type 1 diabetes is the result of insufficient insulin secretion into the bloodstream, most often due to an autoimmune attack on the pancreas glands. Type 2 diabetes is caused by the inability of the surface ligands to adsorb the insulin from the bloodstream. The conventional medicines for diabetes mellitus include sulfonylureas, biguanide, thiazolidinediones, alpha-glucosidase inhibitors, and meglitinide. By February 2022, Severe acute respiratory syndrome coronavirus 2 (SARSCoV2) had infected more than 391 million people worldwide, claiming 5.7 million lives and imposing heavy costs on the healthcare system. The present study aimed to assess the potential use of this non-structural SARS-CoV-2 protein in the treatment of type 2 diabetes mellitus. The nsp10 was structurally aligned with GoDrugBank therapeutic agents, and lixisenatide was found to have the most similar chemical structure. This drug is a glucagon-like peptide-1 (GLP1) receptor agonist used for the treatment of type 2 diabetes mellitus. The best molecular docking energy score for these two proteins was -301.47, and the ligand root mean square deviation was calculated to be 107.93 Å. The molecular dynamics for the stability of the nsp10 and GLP1R binding in triplicate for 150 ns demonstrated that the nsp10-GLP1R remained bound for more than 80 ns. This study indicated that the nsp10 protein can be further studied to be used as an antidiabetic medication.

Integrative Multiomics and Regulatory Network Analyses Uncovers the Role of OAS3, TRAFD1, miR-222-3p, and miR-125b-5p in Hepatitis E Virus Infection.

The hepatitis E virus (HEV) is a long-ignored virus that has spread globally with time. It ranked 6th among the top risk-ranking viruses with high zoonotic spillover potential; thus, considering its viral threats is a pressing priority. The molecular pathophysiology of HEV infection or the underlying cause is limited. Therefore, we incorporated an unbiased, systematic methodology to get insights into the biological heterogeneity associated with the HEV. Our study fetched 93 and 2016 differentially expressed genes (DEGs) from chronic HEV (CHEV) infection in kidney-transplant patients, followed by hub module selection from a weighted gene co-expression network (WGCN). Most of the hub genes identified in this study were associated with interferon (IFN) signaling pathways. Amongst the genes induced by IFNs, the 2'-5'-oligoadenylate synthase 3 (OAS3) protein was upregulated. Protein-protein interaction (PPI) modular, functional enrichment, and feed-forward loop (FFL) analyses led to the identification of two key miRNAs, i.e., miR-222-3p and miR-125b-5p, which showed a strong association with the OAS3 gene and TRAF-type zinc finger domain containing 1 (TRAFD1) transcription factor (TF) based on essential centrality measures. Further experimental studies are required to substantiate the significance of these FFL-associated genes and miRNAs with their respective functions in CHEV. To our knowledge, it is the first time that miR-222-3p has been described as a reference miRNA for use in CHEV sample analyses. In conclusion, our study has enlightened a few budding targets of HEV, which might help us understand the cellular and molecular pathways dysregulated in HEV through various factors. Thus, providing a novel insight into its pathophysiology and progression dynamics.

In silico study of SARS-CoV-2 spike protein RBD and human ACE-2 affinity dynamics across variants and Omicron subvariants.

The coronavirus disease 2019 virus outbreak continues worldwide, with many variants emerging, some of which are considered variants of concern (VOCs). The WHO designated Omicron as a VOC and assigned it under variant B.1.1.529. Here, we used computational studies to examine the VOCs, including Omicron subvariants, and one variant of interest.  Here we found that the binding affinity of human receptor angiotensin-converting enzyme 2 (hACE2) and receptor-binding domain (RBDs) increased in the order of wild type (Wuhan-strain) < Beta < Alpha < OmicronBA.5 < Gamma < Delta < Omicron BA.2.75 < BA.1 < BA.3 < BA.2. Interactions between docked complexes revealed that the RBD residue positions like 452, 478, 493, 498, 501, and 505 are crucial in creating strong interactions with hACE2. Omicron BA.2 shows the highest binding capacity to the hACE2 receptor among all the mutant complexes. The BA.5's L452R, F486V, and T478K mutation significantly impact the interaction network in the BA.5 RBD-hACE2 interface. Here for the first time, we report the His505, an active residue on the RBD forming a salt bridge in the BA.2, leading to increased mutation stability. When the active RBD residues are mutated, binding affinity and intermolecular interactions increase across all mutant complexes. By examining the differences in different variants, this study may provide a solid foundation for structure-based drug design for newly emerging variants.

Omicron SARS-CoV-2 Spike-1 Protein's Decreased Binding Affinity to α7nAChr: Implications for Autonomic Dysregulation of the Parasympathetic Nervous System and the Cholinergic Anti-Inflammatory Pathway—An In Silico Analysis

Omicron is the dominant strain of COVID-19 in the United States and worldwide. Although this variant is highly transmissible and may evade natural immunity, vaccines, and therapeutic antibodies, preclinical results in animal models and clinical data in humans suggest omicron causes a less severe form of infection. The molecular basis for the attenuation of virulence when compared to previous variants is currently not well understood. Using protein–ligand docking simulations to evaluate and compare the capacity of SARS-CoV-2 spike-1 proteins with the different COVID-19 variants to bind to the human α7nAChr (i.e., the core receptor under the control of the vagus nerve regulating the parasympathetic nervous system and the cholinergic anti-inflammatory pathway), we found that 10 out of the 14 mutated residues on the RBD of the B.1.1.529 (Omicron) spike, compared to between 0 and 2 in all previous variants, were present at the interaction interface of the α7nAChr. We also demonstrated, through protein–protein docking simulations, that these genetic alterations cause a dramatic decrease in the ability of the Omicron SARS-CoV-2 spike-1 protein to bind to the α7nAChr. These results suggest, for the first time, that the attenuated nature of Omicron infection in humans and animals compared to previous variants may be attributable to a particular set of genetic alterations, specifically affecting the binding site of the SARS-CoV-2 spike-1 protein to the α7nAChr. [ FROM AUTHOR]

Exploring the binding capacity of lactic acid bacteria derived bacteriocins against RBD of SARS-CoV-2 Omicron variant by molecular simulations.

The changes in the SARS-CoV-2 genome have resulted in the emergence of new variants. Some of the variants have been classified as variants of concern (VOC). These strains have higher transmission rate and improved fitness. One of the prevalent were the Omicron variant. Unlike previous VOCs, the Omicron possesses fifteen mutations on the spike protein's receptor binding domain (RBD). The modifications of spike protein's key amino acid residues facilitate the virus' binding capability against ACE2, resulting in an increase in the infectiousness of Omicron variant. Consequently, investigating the prevention and treatment of the Omicron variant is crucial. In the present study, we aim to explore the binding capacity of twenty-two bacteriocins derived from Lactic Acid Bacteria (LAB) against the Omicron variant by using protein-peptidedocking and molecular dynamics (MD) simulations. The Omicron variant RBD was prepared by introducing fifteen mutations using PyMol. The protein-peptide complexes were obtained using HADDOCK v2.4 docking webserver. Top scoring complexes obtained from HADDOCK webserver were retrieved and submitted to the PRODIGY server for the prediction of binding energies. RBD-bacteriocin complexes were subjected to MD simulations. We discovered promising peptide-based therapeutic candidates for the inhibition of Omicron variant for example Salivaricin B, Pediocin PA 1, Plantaricin W, Lactococcin mmfii and Enterocin A. The lead bacteriocins, except Enterocin A, are biosynthesized by food-grade lactic acid bacteria. Our study puts forth a preliminary information regarding potential utilization of food-grade LAB-derived bacteriocins, particularly Salivaricin B and Pediocin PA 1, for Covid-19 treatment and prophylaxis.Communicated by Ramaswamy H. Sarma.

Human TMPRSS2 non-catalytic ectodomain and SARS-CoV-2 S2' subunit interaction mediated SARS-CoV-2 endocytosis: a model proposal with virtual screening for potential drug molecules to inhibit this interaction.

This study proposes a novel model for integration of SARS-CoV-2 into host cell via endocytosis as a possible alternative to the prevailing direct fusion model. It is known that the SARS-CoV-2 spike protein undergoes proteolytic cleavage at S1-S2 cleavage site and the cleaved S2 domain is primed by the activated serine protease domain (SPD) of humanTMPRSS2 to become S2'. The activated SPD of TMPRSS2 is formed after it is cleaved by autocatalysis from the membrane bound non-catalytic ectodomain (hNECD) comprising of LDLRA CLASS-I repeat and a SRCR domain. It is known that the SRCR domains as well as LDLRA repeat harboring proteins mediate endocytosis of viruses and certain ligands. Based on this, we put forward a hypothesis that the exposed hNECD binds to the S2' as both are at an interaction proximity soon after S2 is processed by the SPD and this interaction may lead to the endocytosis of virus. Based on this hypothesis we have modelled the hNECD structure, followed by docking studies with the known 3D structure of S2'. The interaction interface of hNECD with S2' was further used for virtual screening of FDA-approved drug molecules and Indian medicinal plant-based compounds. We also mapped the known mutations of concern and mutations of interest on interaction interface of S2' and found that none of the known mutations map onto the interaction interface. This indicates that targeting the interaction between the hNECD of TMPRSS2 and S2' may serve as an attractive therapeutic target.Communicated by Ramaswamy H. Sarma.

Receptor binding domain of SARS-CoV-2 from Wuhan strain to Omicron B.1.1.529 attributes increased affinity to variable structures of human ACE2.

BACKGROUND: COVID-19 is an infectious disease declared as a global pandemic caused by SARS-CoV-2 virus. Genomic changes in the receptor binding domain (RBD) region of SARS-CoV-2 led to an increased, infectivity in humans through interaction with the angiotensin-converting enzyme2 (ACE2) receptor. Simultaneously, the genetic variants in ACE2 provide an opportunity for SARS-CoV-2 infection and severity. We demonstrate the binding efficiencies of RBDs of SARS-CoV-2 strain with ACE2 variants of the human host. METHODOLOGY: A Total of 615 SARS-CoV-2 genomes were retrieved from repository. Eighteen variations were identified contributing to structural changes in RBD that are distributed in 615 isolates. An analyses of 285 single nucleotide variances at the coding region of the ACE2 receptor showed 34 to be pathogenic. Homology models of 34 ACE2 and 18 RBD structures were constructed with 34 and 18 structural variants, respectively. Protein docking of 612 (34 *18) ACE2-RBD complexes showed variable affinities compared to wildtype Wuhan's and other SARS-CoV-2 RBDs, including Omicron B.1.1.529. Finally, molecular dynamic simulation was performed to determine the stability of the complexes. RESULTS: Among 612, the top 3 complexes showing least binding energy were selected. The ACE2 with rs961360700 variant showed the least binding energy (-895.2 Kcal/mol) on binding with the RBD of Phe160Ser variant compared to Wuhan's RBD complex. Interestingly, the binding energy of RBD of Omicron B.1.1.529 with ACE2 (rs961360700) structure showed least binding energy of -1010 Kcal/mol. Additionally, molecular dynamics showed structure stability for all the analysed complexes with the RMSD (0.22-0.26 nm), RMSF (0.11-0.13 nm), and Rg (2.53-2.56 nm). CONCLUSION: In conclusion, our investigation highlights the clinical variants contributing to structural variants in ACE2 receptors that lead to efficient binding of SARS-CoV-2. Therefore, screening of these ACE2 polymorphisms will help detect COVID-19 risk population so as to provide additional care and for safe management.

Virtual screening by targeting proteolytic sites of furin and TMPRSS2 to propose potential compounds obstructing the entry of SARS-CoV-2 virus into human host cells.

BACKGROUND AND AIM: The year 2020 begins with the outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that cause the disease COVID-19, and continue till today. As of March 23, 2021, the outbreak has infected 124,313,054 worldwide with a total death of 2,735,707. The use of traditional medicines as an adjuvant therapy with western drugs can lower the fatality rate due to the COVID-19. Therefore, in silico molecular docking study was performed to search potential phytochemicals and drugs that can block the entry of SARS-CoV-2 into host cells by inhibiting the proteolytic cleavage activity of furin and TMPRSS2. EXPERIMENTAL PROCEDURE: The protein-protein docking of the host proteases furin and TMPRSS2 was carried out with the virus spike (S) protein to examine the conformational details and residues involved in the complex formation. Subsequently, a library of 163 ligands containing phytochemicals and drugs was virtually screened to propose potential hits that can inhibit the proteolytic cleavage activity of furin and TMPRSS2. RESULTS AND CONCLUSION: The phytochemicals like limonin, gedunin, eribulin, pedunculagin, limonin glycoside and betunilic acid bind at the active site of both furin and TMPRSS2. Limonin and gedunin found mainly in the citrus fruits and neem showed the highest binding energy at the active site of furin and TMPRSS2, respectively. The polyphenols found in green tea can also be useful in suppressing the furin activity. Among the drugs, the drug nafamostat may be more beneficial than the camostat in suppressing the activity of TMPRSS2.

Investigation of the binding and dynamic features of A.30 variant revealed higher binding of RBD for hACE2 and escapes the neutralizing antibody: A molecular simulation approach.

With the emergence of Delta and Omicron variants, many other important variants of SARS-CoV-2, which cause Coronavirus disease-2019, including A.30, are reported to increase the concern created by the global pandemic. The A.30 variant, reported in Tanzania and other countries, harbors spike gene mutations that help this strain to bind more robustly and to escape neutralizing antibodies. The present study uses molecular modelling and simulation-based approaches to investigate the key features of this strain that result in greater infectivity. The protein-protein docking results for the spike protein demonstrated that additional interactions, particularly two salt-bridges formed by the mutated residue Lys484, increase binding affinity, while the loss of key residues at the N terminal domain (NTD) result in a change to binding conformation with monoclonal antibodies, thus escaping their neutralizing effects. Moreover, we deeply studied the atomic features of these binding complexes through molecular simulation, which revealed differential dynamics when compared to wild type. Analysis of the binding free energy using MM/GBSA revealed that the total binding free energy (TBE) for the wild type receptor-binding domain (RBD) complex was -58.25 kcal/mol in contrast to the A.30 RBD complex, which reported -65.59 kcal/mol. The higher TBE for the A.30 RBD complex signifies a more robust interaction between A.30 variant RBD with ACE2 than the wild type, allowing the variant to bind and spread more promptly. The BFE for the wild type NTD complex was calculated to be -65.76 kcal/mol, while the A.30 NTD complex was estimated to be -49.35 kcal/mol. This shows the impact of the reported substitutions and deletions in the NTD of A.30 variant, which consequently reduce the binding of mAb, allowing it to evade the immune response of the host. The reported results will aid the development of cross-protective drugs against SARS-CoV-2 and its variants.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Molecular Docking and Dynamics Studies to Explore Effective Inhibitory Peptides Against the Spike Receptor Binding Domain of SARS-CoV-2.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a pandemic due to the high transmission and mortality rate of this virus. The world health and economic sectors have been severely affected by this deadly virus, exacerbated by the lack of sufficient efficient vaccines. The design of effective drug candidates and their rapid development is necessary to combat this virus. In this study, we selected 23 antimicrobial peptides from the literature and predicted their structure using PEP-FOLD 3.5. In addition, we docked them to the SARS-CoV-2 spike protein receptor-binding domain (RBD) to study their capability to inhibit the RBD, which plays a significant role in virus binding, fusion and entry into the host cell. We used several docking programs including HDOCK, HPEPDOCK, ClusPro, and HawkDock to calculate the binding energy of the protein-peptide complexes. We identified four peptides with high binding free energy and docking scores. The docking results were further verified by molecular dynamics (MD) simulations to characterize the protein-peptide complexes in terms of their root-mean-square fluctuation (RMSF), root-mean-square deviation (RMSD), radius of gyration (Rg), solvent-accessible surface area (SASA), and hydrogen bond formation. Allergenicity and toxicity predictions suggested that the peptides we identified were non-allergenic and non-toxic. This study suggests that these four antimicrobial peptides could inhibit the RBD of SARS-CoV-2. Future in vitro and in vivo studies are necessary to confirm this.

Structural Analysis on the Severe Acute Respiratory Syndrome Coronavirus 2 Non-structural Protein 13 Mutants Revealed Altered Bonding Network With TANK Binding Kinase 1 to Evade Host Immune System.

Mutations in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have made this virus more infectious. Previous studies have confirmed that non-structural protein 13 (NSP13) plays an important role in immune evasion by physically interacting with TANK binding kinase 1 (TBK1) to inhibit IFNß production. Mutations have been reported in NSP13; hence, in the current study, biophysical and structural modeling methodologies were adapted to dissect the influence of major mutations in NSP13, i.e., P77L, Q88H, D260Y, E341D, and M429I, on its binding to the TBK1 and to escape the human immune system. The results revealed that these mutations significantly affected the binding of NSP13 and TBK1 by altering the hydrogen bonding network and dynamic structural features. The stability, flexibility, and compactness of these mutants displayed different dynamic features, which are the basis for immune evasion. Moreover, the binding was further validated using the MM/GBSA approach, revealing that these mutations have higher binding energies than the wild-type (WT) NSP13 protein. These findings thus justify the basis of stronger interactions and evasion for these NSP13 mutants. In conclusion, the current findings explored the key features of the NSP13 WT and its mutant complexes, which can be used to design structure-based inhibitors against the SARS-CoV-2 new variants to rescue the host immune system.

GRP78: A possible relationship of COVID-19 and the mucormycosis; in silico perspective.

Mucormycosis is a severe fungal infection reported in many cancer survivors, diabetic and immune-suppressed patients during organ transplants. A vast spark in the reported COVID-19 cases is noticed in India during the second wave in May 2021, when Mucormycosis is declared an epidemic. Despite being a rare disease, the mortality rate associated with Mucormycosis is more than 40%. Spore coat proteins (CotH) are essential proteins in many pathogenic bacteria and fungi. CotH3 was reported as the vital protein required for fungal virulence in Mucormycosis. We previously reported the involvement of the host cell-surface receptor GRP78 in SARS-CoV-2 spike recognition. Additionally, GRP78 is known to be the virulence factor during Mucormycosis. Using state-of-the-art structural bioinformatics and molecular modeling tools, we predicted the GRP78 binding site to the Rhizopus delemar CotH3 protein. Our findings pave the way toward rationally designing small molecule inhibitors targeting the GRP78 and its counter proteins in both pathogenic viral (SARS-CoV-2 spike) and fungal (R. delemar CotH3) diseases.

Peptide Platform as a Powerful Tool in the Fight against COVID-19.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in a global pandemic causing over 195 million infections and more than 4 million fatalities as of July 2021.To date, it has been demonstrated that a number of mutations in the spike glycoprotein (S protein) of SARS-CoV-2 variants of concern abrogate or reduce the neutralization potency of several therapeutic antibodies and vaccine-elicited antibodies. Therefore, the development of additional vaccine platforms with improved supply and logistic profile remains a pressing need. In this work, we have validated the applicability of a peptide-based strategy focused on a preventive as well as a therapeutic purpose. On the basis of the involvement of the dipeptidyl peptidase 4 (DPP4), in addition to the angiotensin converting enzyme 2 (ACE2) receptor in the mechanism of virus entry, we analyzed peptides bearing DPP4 sequences by protein-protein docking and assessed their ability to block pseudovirus infection in vitro. In parallel, we have selected and synthetized peptide sequences located within the highly conserved receptor-binding domain (RBD) of the S protein, and we found that RBD-based vaccines could better promote elicitation of high titers of neutralizing antibodies specific against the regions of interest, as confirmed by immunoinformatic methodologies and in vivo studies. These findings unveil a key antigenic site targeted by broadly neutralizing antibodies and pave the way to the design of pan-coronavirus vaccines.

Potential inhibitors for the novel coronavirus (SARS-CoV-2).

The lack of a vaccine or any effective treatment for the aggressive novel coronavirus disease (COVID-19) has created a sense of urgency for the discovery of effective drugs. Several repurposing pharmaceutical candidates have been reported or envisaged to inhibit the emerging infections of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), but their binding sites, binding affinities and inhibitory mechanisms are still unavailable. In this study, we use the ligand-protein docking program and molecular dynamic simulation to ab initio investigate the binding mechanism and inhibitory ability of seven clinically approved drugs (Chloroquine, Hydroxychloroquine, Remdesivir, Ritonavir, Beclabuvir, Indinavir and Favipiravir) and a recently designed α-ketoamide inhibitor (13b) at the molecular level. The results suggest that Chloroquine has the strongest binding affinity with 3CL hydrolase (Mpro) among clinically approved drugs, indicating its effective inhibitory ability for SARS-CoV-2. However, the newly designed inhibitor 13b shows potentially improved inhibition efficiency with larger binding energy compared with Chloroquine. We further calculate the important binding site residues at the active site and demonstrate that the MET 165 and HIE 163 contribute the most for 13b, while the MET 165 and GLN 189 for Chloroquine, based on residual energy decomposition analysis. The proposed work offers a higher research priority for 13b to treat the infection of SARS-CoV-2 and provides theoretical basis for further design of effective drug molecules with stronger inhibition.